clc;
syms u v;
syms a positive;
x=u*cos(v);
y=u*sin(v);
z=v;

f=x^2*y+z*y^2;

I=surf_integral(f,[x,y,z],[u,v],[0,a],[0,2*pi]);
pretty(I)
latex(I)